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Abstract. We analyse the respective benefits and drawbacks of ground-based and space-based transit surveys 
for extrasolar planets. Based on simple but realistic assumptions about the fraction of lower main sequence stars 
harboring telluric and giant planets within the outer limit of the habitable zone, we predict the harvests of fictitious 
surveys with three existing wide field optical and near-IR cameras: the CFHT-Megacam, SUBARU-Suprime and 
VISTA-IR. An additional promising instrument is considered, VISTA-Vis, currently under development. The 
results are compared with the harvests predicted under exactly the same assumptions, for the space missions 
COROT and KEPLER. We show that ground-based wide field surveys may discover more giant planets than 
space missions. However, space surveys seem to constitute the best strategy to search for telluric planets. In this 
respect, the KEPLER mission appears 50 times more efficient than any of the ground-based surveys considered 
here. KEPLER might even discover telluric planets in the habitable zone of their host star. 
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1. Introduction 

More than 160 extrasolar planets have been discovered^ 
since the first detection of a "Hot Jupiter" around the 
main-sequence star 51 Peg | |Mayor fc Queloz, 1995| ). Most 
of these discoveries are the results of radial velocity (RV) 
searches. However, while RV searches are very efficient 
at detecting massive planets, they do not allow, on their 
own, to characterize completely the orbit of the planet. 
Due to the degeneracy between the mass of the planet 
and the inclination of its orbit with respect to the line 
of sight (see, |e.g., Perryman, 2000, for a review] ), M smi 
is the only measurable quantity using the RV technique 
alone, where AI is the mass of the planet and i, the incli- 
nation of its orbit. At the end of 1999, a significant advance 
was accomplished by the observation of the first transit of 
an extrasolar planet discovered through the RV method 
HHenry et al., 2000| |Charbonneau et aL, 2000| ). With si- 
multaneous spectroscopic and photometric observations of 
the star during the transit, both the mass and the radius 
of the companion object could be determined, confirming 
beautifully its planetary nature, and demonstrating the 
great interest of combining RV and transit observations. 
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Transit photometry is not only of interest as a fol- 
low up of RV searches, but also as a way to find 
new planets. It may prove even more sensitive to low 
mass planets than RVs. In fact, the method is con- 
sidered as one of the most promising ways of finding 
earth-class planets ( |Sackett, 1999) [Schneider, 2000| ). The 
OGLE survey (Optical Gravitational Lensing Experiment, 
[Udalski et al., 2002al 12002 b|l , a ground-based experiment 
initially devoted to stellar microlensing in the galac- 
tic bulge, has been extended to the search for plane- 
tary transits. In 2003, Konacki et al. l|7003a, 2003b1 an- 
nounced the characterization of the first exoplanet discov- 
ered by transit photometry, orbiting the star OGLE-TR- 
56. The surprisingly short period of 1.2 day of the com- 
panion object was confirmed by RV follow-up, a value 
much below the lower end of the period distribution 
of planets detected by RV surveys ( |Udry et al., 2003| ). 
This important discovery demonstrated the effective- 
ness of transit photometry at detecting extrasolar plan- 
ets, and was soon followed by four other discoveries 
by the same team ||Bouchy et al., 2004||Pont et al., 2004| 
[Konacki et al., 2005| |. Other surveys have started, such as 
TrES (Trans-Atlantic Exoplanet Survey) multisite tran- 
sit survey, and its first discovery of an exoplanet, TrES- 
1 | |Alonso et al., 200^ . Very recently, transits of a satur- 
nian size planet discovered through the RV method have 
been observed by the N2K Consortium ( [Sato et al., 2005] ) 
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around the star HD 149026, leading to a total of eight 
stars known to experience planetary transits. 

With the potential of the transit method now well 
established, space missions have been proposed in the 
hope to discover dozens of new planets. Among them 
are the COROT ( |Rouan et al., 2000| ) and KEPLER 
dKoch et al., 1998| ) missions, designed to perform high 
precision photometry of selected fields over several years. 
In parallel, wide field cameras are becoming available on 
large ground-based telescopes, making it possible to sur- 
vey, from the ground, large numbers of stars simultane- 
ously. The aim of the present paper is to weight the rel- 
ative benefits and drawbacks of ground-based and space 
surveys, by carrying out a comparative study of their effi- 
ciency based on simple simulations. In particular, we in- 
vestigate whether a large dedicated ground-based survey 
could be a competitive alternative to a high-cost space 
mission, or whether both approaches are complementary. 

On the basis of realistic assumptions, we estimate the 
harvest in extrasolar planets for each type of survey. Our 
main goal is to compare the merits of each survey type, 
rather than predicting absolute discovery rates. Such ab- 
solute predictions are far too sensitive to the assumptions 
about the physical properties and about the formation 
process of planetary systems. We consider only simple dis- 
tributions of planets around main sequence stars, in terms 
of radii and orbital distances. These are the only relevant 
quantities needed to carry out a comparative study, along 
with estimates of the density of stars in the simulated 
stellar fields. 

The method used to estimate the harvests in planets 
is described in Sect. 2. The surveys considered and the 
results of our simulations are presented in Sect. 3, while 
the results and conclusions are summarized in Sect. 4. 

2. Description of the method 

2.1. Strategy and field selection 

Planetary transits are rare and the eclipse seen in their 
host stars is extremely shallow: around 10"^ mag for a 
Jupiter-size planet and 10"'* mag for an Earth-size planet 
in orbit around of a solar-type star. The question of the 
compromise that inevitably has to be made between field 
and depth is then a critical issue for any transit search. 
In addition, the time scales involved, i.e. the (expected) 
duration of the transits, is limiting the range of suitable 
exposure times, as soon as the goal is to properly sample 
the light curves during the transit phase. The maximum 
exposure for a given telescope is therefore mainly dictated 
by the size of the planets one wished to discover with re- 
spect to the size of its host star. Clearly, large telescopes 
are required in order to apply as short as possible expo- 
sure times and to sample the light curves at best. Wide 
field is also of major importance to observe many stars 
simultaneously and to compensate for the scarcity of tran- 
sit events (e.g., |Mallen-Ornelas et al., 2003| ). Being aware 
of this general line to design surveys, we describe in the 



following the simulations used in order to estimate the 
harvests of fictitious and real planetary transit surveys. 

2.2. Suitable host stars 

2.2.1. Density and spectral types 

For a given planetary radius, the transit depth is inversely 
proportional to the square of the stellar radius. Thus, plan- 
etary transits will be more easily detected in cool dwarf 
stars. In this paper, we only consider the main-sequence 
stars with spectral subtypes from FO to A/9 as suitable 
candidates. We refer to them as Lower Main-Sequence 
Stars (LMSS). 

The galactic plane is the obvious place to look at, in or- 
der to observe simultaneously many LMSS. It is, indeed, 
the strategy adopted in the pioneering work by Mallen- 
Ornelas et al. (2003). Since the highest possible Signal-to- 
Noise Ratio {SNR) per exposure is required, fields with 
minimal reddening are prefered. A best choice consists of 
fields close to, but not right in, the galactic plane, with 
galactic latitudes between 2° and 6°. We have chosen to 
carry out all our computations for fields with a mean lati- 
tude of 4° , where a representative extinction coefficient is 
Av=0.7 mag/kpc ( .Schlegel et al., 1998| ). 

The projected densities of LMSS at low galactic lat- 
itude are estimated from the Gliese & Jarheiss II1991|I 
and the Zakhozhaj H1979|l catalogues of stars in the solar 
neighborhood (see Fig. We assume that these stellar 
densities also apply to the rest of the galactic plane. This 
is probably a fair approximation as long as the galactic 
bulge is avoided. The fields we are modeling here mainly 
contain stars belonging to the spiral arms of our Galaxy. 

Under these assumptions, the number of stars in the 
field of view is computed for every half-magnitude bin 
and for all spectral subtypes from FO to M9. The vol- 
umes sampled by the Gliese et al. (1991) and Zakhozhaj 
et al. (1979) catalogues are scaled to the volumes of each 
of the surveys considered. Thus, the total number of stars 
observed per half-magnitude bin depends on the size of 
the field of view and the filter used. Fig. |21 shows the dis- 
tribution of LMSS, computed for all subspectral types 
and for 4 limiting magnitudes in the /-band. Going to 
the infrared allows to observe a larger fraction of very- 
low-mass stars, as is shown in Fig. |21 

2.2.2. Blends and binarity 

All our calculations are exclusively made for LMSS and 
do not include the effect of blends, for two reasons: 

1. Introducing all spectral types will increase the num- 
ber of potentially interesting transits. It will also 
increase the crowding of the fields. Since stars other 
than LMSS have larger diameters, the net result of 
the competition between the two effects is likely a 
reduction of the true observable number of transits. 
Our results are therefore upper limits on the number 
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FO F5 GO G5 KO K5 MO M5 

Spectral subtype 



Fig. 1. Density of Lower Main Sequence Stars (LAISS) 
spectral subtypes in the solar neighborhood, estimated 
from the Gliese (1991) and Zakhozhaj (1979) nearby stars 
catalogues. 
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Fig. 2. Spectral subtypes distribution of LAISS for 4 
limiting magnitudes (15, 17, 19 and 21). An extinction 
coefficient Av=0.7 mag/kpc and a maximum distance 
Dl=4300 pc are used (see text). 

of transits seen by each survey. However, since our 
goal is a comparative study between surveys, rather 
than an attempt to predict absolute transit counts. 



FO F5 GO G5 KO K5 MO M5 FO F5 GO G5 KO K5 MO M5 

Spectral subtype Spectral subtype 

Fig. 3. Same Figure as Fig. |21 but in the if-band. The 
proportion of very-low-mass stars, for a given limiting 
magnitude, is larger than in the /-band. 

leaving aside the spectral types other than FO to AI9 
does not affect our conclusions. 

2. More blends also mean an increased number of false de- 
tections. For Jovian-mass planets, a large part of these 
false detections can be rejected through a fine-tuned 
transit light-curve analysis or with follow-up RV ob- 
servations dBouchy et al., 2005| ). In some cases, decid- 
ing whether an eclipse is due to a genuine planet or not 
can be very tricky (see e.g. [Mandushev et al., 2005| ). 
The difficulty increases drastically in the case of a tel- 
luric planet. It does not only depend on the planet 
type, but also on the method of analysis used to build 
the stellar light curves and on the transit detection 
method and selection criteria: methods based on aper- 
ture photometry will be heavily affected by seeing and 
blends, while differential imaging, PSF fitting, or im- 
age deconvolution behave very differently with respect 
to blends and to PSF mismatch. Finally, prior knowl- 
edge on the exact field geometry can be introduced 
as well in the analysis, e.g., by targetting fields with 
archived HST images used as reference maps to carry 
out the photometry. The degree of additional compli- 
cation due to blends is such that it falls well out of 
the scope of the present paper. Again, only relative 
estimates between surveys are aimed at. 

Very little is known about the formation of planets in 
binary systems. Furthermore, detecting transits in close 
binary systems is possible but can be difficult, because 
of the decrease of the transit depth. We have chosen to 
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adopt an arbitrary fraction of 50% of binary LMSS and 
to do not allow a transit detection around them. The exact 
value of this fraction has no influence on the result of our 
comparative study. Our whole simulations are designed to 
predict the maximum number of planets that can poten- 
tially be discovered with each type of survey, with a slight 
and arbitrary correction for binarity. 

2.2.3. Distance 

Among our assumptions is the maximum distance above 
which the density of LMSS may not be valid. As we con- 
sider fields with a low but non-zero galactic latitude, this 
maximum distance can be computed from the thickness of 
the galactic disk. Considering a thin disk with a thickness 
of 600 pc, positioning the Sun at mid height, and using 
an average galactic latitude of 4° for all the surveys con- 
sidered in this work, the limiting distance Dl is simply: 

In the following we will consider that the density of 
LMSS drops to zero beyond the distance D^. The limit- 
ing magnitudes imposed to us by the telescopes and cam- 
eras chosen in this work for ground-based surveys are 
well above the limiting magnitude due to pure distance 
effects. The magnitude range is therefore limited by dis- 
tance rather than by instrumental depth. For space-based 
surveys, the instrumental depth is the limiting factor. 
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Fig. 4. Top: Radii distribution of giant planets interpo- 
lated from the four giant planets in our Solar System. 
The y-axis gives the probability that a planet has a ra- 
dius larger than the value on the x-axis, in earth radii. 
Bottom: Radii distribution of THZ (solid line) and TIZ 
(dotted line) , interpolated from the simulations of the for- 
mation of telluric planets by Chambers (2001). 



2.3. Estimated number of planets in the field 

Once the number of target stars has been estimated, one 
also needs to make assumptions about the average number 
of planets hosted by LMSS, and on the distribution of 
their radii and periods. We consider six types of planets: 

— Very Hot Jupiters (VHJ): we define these objects 
as giant planets with a maximum period of 3 days and 
an arbitrary minimum period of 1 day. On the basis 
of the results of Gaudi et al. (2005)] we assume that 
0.2% of LMSS may host a VHJ. We use a rough 
radii distribution represented by a sigmoid function, 
interpolated from the radii distribution observed in 
our own solar system (see Fig. 0}. We define the 
V H J-Zone as the range of distances to the host star 
corresponding to orbital periods between 1 and 3 
days. We assume that the periods of the VHJ are 
uniformly distributed in the VHJ-Zone. 

— Hot Jupiters (HJ): are defined as giant planets 
with a period in between 3 and 9 days. Following 
Gaudi et al. (2005), we adopt a value of 1% for the 
fraction of LMSS orbited by such objects. We use 
the same radii distribution as for the VHJ. We define 
the HJ-Zone as the region around a star where the 
period of a planet is between 3 and 9 days. As noticed 



by Gaudi et al. (2005), about half of the HJ have a 
period lower than 3.5 days. We thus consider that half 
of the HJ have periods uniformly distributed in the 
range 3-3.5 days, and that the other half have periods 
uniformly distributed in the range 3.5-9 days. HJ 
and VHJ forms the group of the Close-in Extrasolar 
Giant Planets {CEGP). We shall denote them {y)HJ 
when considered as a single group. 

Giant planets in the Habitable Zone {GHZ): 
giant planets in the zone identified by Kasting et al. 
(1993) as suitable for life. We use the limits estab- 
lished by Kasting et al. (1993) for the Zero Age Main 
Sequence HZ {ZAMS HZ) and their "intermediate" 
habitability criteria. Based on the measured fraction 
of extrasolar planets found in the HZ, we assume 
that 3% of the LMSS harbor such a planet. We use 
the same radii distribution for these objects as for the 
HJ and VHJ. We consider that the periods of the 
GHZ are uniformly distributed across the HZ . 

Giant planets in the Intermediate Zone (GIZ): 

the Intermediate Zone is located between the HZ and 
the HJ-Zone. We assume that the fraction of LMSS 
with a GIZ is the same as in the HZ, i.e., 3%. The 
same radii distribution is used for GIZ as for the 
other types of giants planets considered. We consider 
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that the periods of the GIZ are distributed uniformly 
in the Intermediate Zone. 

- Telluric planets in the Habitable Zone (THZ): 
nothing is known about telluric planets outside our 
own Solar System. We base our estimates on the 
terrestrial planets formation simulations of Chambers 
{2001). Using these results and the assumption that 
all LM SS are surrounded by a protoplanetary disk 
at the beginning of their life (which will lead to the 
formation of telluric planets) we obtain 1.78 as the 
average number of terrestrial planets in the HZ per 
LMSS. This number does not take into account the 
fact that, for a poorly constrained fraction of the stars, 
one or more giant planet may migrate inwards over 
the outer limit of the HZ , hence affecting the orbits 
of the telluric planets and maybe accreting them or 
expulsing them on a hyperbolic orbit. We take this 
fraction as the sum of the fractions of VHJ, HJ, 
GHZ and GIZ, leading to 7% of the LMSS hosting 
a giant planet inside the outer limit of the HZ. 
Assuming that no terrestrial planet can survive under 
these conditions, we obtain 1.78 x (1 — 0.07) « 1.66 
as the average number of terrestrial planets in the 
HZ per LMSS. We also use the results of Chambers 
H2()()l|l . interpolated with a sigmoid function (see 
Fig. ^ to estimate the radii distribution of telluric 
planets in the HZ. We consider that the periods of 
the THZ are distributed uniformly in the HZ . 

— Telluric planets in the Intermediate Zone 

(TIZ): using the results of Chambers H2()()l|l and 
correcting for the fraction of LMSS having a giant 
planet inside the outer limit of the IZ (4%), we 
estimate the average number of terrestrial planets in 
the IZ to be 1.51 per LMSS. The results of Chambers 
are again used to estimate the radii distribution of 
telluric planets in the IZ (Fig. 0)). We consider that 
the periods of the TIZ are distributed uniformly in 
the IZ. 



The planets outside the HZ are not considered here, 
as their long period and their low geometric probability 
of producing a detectable transit make them very poor 
candidates. We also do not consider telluric planets with a 
period shorter than 9 days, because they are not predicted 
by the simulations of Chambers and because there is no 
such planet in our own solar system. 

Our assumptions about giant planets are empirical, 
and based on previous surveys that are certainly biased 
towards a given type of planet and orbit and still suffer 
from low number statistics. Furthermore, we assume that 
the average expected number of planets for the four zones 
we defined earlier are the same for all spectral subtypes 
from FO to M9. This assumption does not have a strong 
influence on the results given our already rough knowl- 
edge of planetary formation, although a metallicity de- 



Table 1. Instruments considered in our fictitious ground- 
based surveys. 



Instr. 


Telescope 


Field 


Location 




diameter 


(sq. degrees) 




Megacam 


3.6 m 


0.90 


Hawaii (CFHT) 


Suprime 


8.3 m 


0.26 


Hawaii (Subaru) 


Vista-IR 


4.0 m 


1.00 


Chile (VLT) 


Vista-Vis 


4.0 m 


2.25 


Chile (VLT) 



pendence on exoplanets population has been emphasized 
(see e.g. [Santos et al., 2004| ) and low-mass stars tend to 
be older on average and have a lower metallicity. The pre- 
dicted harvests should nevertheless be of the right order 
of magnitude and are certainly adequate for comparison 
purposes. 

The assumptions made on telluric planets are based 
on a purely theoretical work. We do not take into account 
the possible existence of giant planet cores formed in the 
outer region of the disk and that later migrated inwards. 
These objects, without a massive primary atmosphere, 
could consist of large planets of pure rock and ice. We 
prefer to remain conservative and to consider exclusively 
planets formed in the more standard scenario thought to 
be responsible for the formation of planetary systems such 
our own. 

The known weaknesses of the assumptions used in our 
simulations are not a critical issue for our purpose, which 
is to compare the relative strengths and weaknesses of 
different transit surveys rather than to predict actual dis- 
covery rates for a given survey. Our results may be scaled 
up and down, but we believe that they remain correct as 
long as the aim is to carry out comparative studies. 

2.4. Parameters of the instruments cliosen for tfie 
fictitious ground-based surveys 

Transit searches must combine large field of view, depth 
and good temporal sampling. Although the range of pos- 
sible telescope/camera combinations is broad, we have se- 
lected four instruments that can be considered as represen- 
tative of the present or soon available astronomical facili- 
ties. Three of these instruments are already in use: CFHT- 
Megacam, Subaru-Suprime, and Vista-IR. The latter one 
is only proposed: Vista- Vis^ (see Table 01. 

For each instrument, we have estimated the SNR and 
the saturation magnitude for a range of exposure times, for 
a fixed airmass (1.6) and typical seeing (1 arcsec). These 
parameters are computed for different filters, using the 
Exposure Time Calculators (ETC) available for Subaru 
and the CFHT. The bright cut applied to the magnitude 
distribution of host stars corresponds to saturation time. 
The faint cut corresponds to the magnitude of stars that 
have a SNR < 10, where a shallow transit would not be 
detected. 

^ see VISTA web site, 
http: / /www. vista.ac.uk( 
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For Vista~Vis, we used the information given in the 
online ETC of the EMMI camera. This ESO instrument 
is mounted on the 3.5 m NTT at La Silla (Chile) and has a 
throughput similar to Vista, which will also be a 4m-class 
telescope. For Vista-IR, we use the information given on 
the online ETC of the SOFI IR camera of the ESO 3.5m 
NTT. 

All our estimated SNR are corrected for the extra 
source of noise introduced by stellar variability, following 
the formula: 



SNR2 



SNRi 



^/l + ia,SNRi) 



(2) 



where SNRi and SNR2 are the SNR before and 
after taking cr* into account, the standard deviation of 
the stellar variability. We adopt cr, = 100 ppm as esti- 
mated from the solar variability within a frequency in- 
terval that matches our adopted range of transit dura- 
tions ( |Borde, 2003| ). We consider this same value for every 
LMSS spectral subtypes. 

2.5. Computation of the harvests and detection criteria 

Under all the above assumptions, the expected numbers 
of transit detections are estimated as follows. 

For every spectral subtype FO to M9, we compute 100 
circular orbits for each of the four zones considered, dis- 
tributed as mentioned in Section 2.3. The geometric tran- 
sit probability for an orbit is given by 



p _ ^* 



(3) 



where i?* is the radius of the star and a the semi- 
major axis of the planet orbit. For every orbit and for a 
range of inclinations and planetary radii, we compute the 
total crossing-time of the transit and the duration of the 
flat part of the light curve, following the calculation by 
Mallen-Ornelas et al. (2001). The results are then aver- 
aged on the inclinations and radii. 

A window function is computed for each survey, based 
on the total number of nights in the campaign and the 
visibility of the field each night, in the case of ground- 
based surveys. To take weather effects into account, we 
assume that for 1 night out of 10, no observations are 
taken at all. In addition, 10 chunks of 1 hour are ran- 
domly removed from the 9 remaining nights to account 
for technical problem, clouds, or unexpected overheads. 

We then compute a probability PvIsN that a transit is 
observed N times for a specific orbit and a specific LMSS. 
We compute this probability for N varying from I to X, 
where X is the maximum number of transits we could ob- 
serve for the same star during the whole observing season. 
As the shortest period considered is 1 day, X is simply the 
duration of the survey in days. 

For a specific spectral type and a specific distance, the 
probability PotsN that a transit occurs N times during the 
survey and is observable is PotsN = Ptr x PvisN- 



We then consider that the observed dimming of a light 
curve can be attributed to a genuine planet only if at least 
three eclipses are detected. We define k as the number of 
transits observed during the survey and impose that the 
SNR of the light curve, integrated on the duration of the 
flat part of the k transits is at least /3 time greater than 
the inverse of transit depth, i.e. 



(4) 



where i?p is the planetary radius and i?* the radius of 
the host star. 

In other words, the significance of planetary transits 
in a stellar light curve is /3 x cr. 

We have chosen to adopt /3 = 9. This value is high 
enough to reject the majority of the statistical artifacts. 
Is is also the one adopted by the OGLE-III team (Udalski, 
private communication). The value /3 = 7 is also tested, 
as it is used in other transit studies such as in Borde et 
al. (2003) for the COROT mission. 

For a specific planet transiting k times during the sur- 
vey, and for a star of a given magnitude and a given radius 
i?,, equation (4) allows a computation of the minimum 
planetary radius needed to get a SNR high enough to al- 
low a statistically significant detection. The relevant radii 
distribution then leads to a determination of the proba- 
bility that the planet has a radius at least equal to the 
minimum radius of detection, i.e. the fraction of planets 
which would produce a dimming strong enough to be de- 
tected. 

For every half magnitude bin and for every zone con- 
sidered, the results are averaged over 100 orbits and multi- 
plied by the number of stars of each spectral type present 
in each half-magnitude bin. We then multiply by the frac- 
tion of stars expected to host a planet (giant or telluric) 
in the zone considered, leading to the expected number of 
planet detections. As a final step, we sum the numbers of 
planets for each half-magnitude bin and spectral subtype, 
and obtain the total number of detections for each planet 
type. 



3. Results 

3.1. Analysis of existing surveys 
EXPLORE-I 



OGLE-III and 



We have tested our simulations 
isting ground-based surveys, 
Mallen-Ornelas et al., 2003D 



on two ex- 
EXPLORE-I 
and OGLE-III 
lUdalski et al., 2002a| I2002bl IMUjl . The goal was 
to compare our predictions to the actual results of these 
surveys and to check the validity of our assumptions. Our 
main goal remains the comparison of existing ground- 
based surveys to future fictitious ground-based surveys 
and to space missions. 
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Table 2. Results predicted for EXPLORE-I. In the first 
line, we adopt Pyis2 =0.1 for {V)HJ and Pvis2 = for the 
other types of planets. The second and third lines involve 
a PvisN computed assuming good weather conditions (see 
text). The values in brackets are for /? = 7, and the others 
correspond to = 9. 



3.2. EXPLORE-I. 

The Extrasolar Planet Occultation Research search lasted 
11 nights, at the CTIO 4m telescope, in the /-band with 
the mosaic II camera. It concentrated on one single 0.36 
deg2 field near the galactic plane {I = -27.8°, b=-2.7°) 
containing - 100,000 stars down to 7=18.2 and ~ 350,000 
stars down to 7=21.0. Poor weather conditions led to a 
degraded window function, affecting the probability P.uis2 
to detect two transits of the same planet. Pyis2 was esti- 
mated as 0.1 for (V^)77 J by the EXPLORE team (Mallen- 
Ornelas, private communication). 

The typical exposure time was 60 seconds, and the de- 
tector read time plus overhead amounted to 101 seconds. 
The analysis of the data is not completed, and no final 
detection criterion has been decided up to now (Mallen- 
Ornelas, private communication). No exoplanet transit 
has been discovered so far. 

Following the simulations presented in this paper, with 
Pvis2 = 0.1 for [V)HJ (and PvisN>2 set to 0, i.e. assum- 
ing that the probability to detect three transits of the 
same planet is negligible), we compute the expected har- 
vest for EXPLORE-I, given its true weather conditions. 
We also compute the harvest for good weather conditions, 
i.e. without fixing PyisN but estimating it according to 
the method described in Section 2.5. Our results for the 
two detection criteria (/3 = 9 and (3 = 1) are presented in 
Table m 

We present in Fig. [S] the magnitude distribution of the 
host stars {V)HJ expected according to our simulation 
under good weather conditions and /? = 9. Fig. [S] also 
shows the distribution of the distances, of the spectral 
subtypes and of the number of transits observed for the 
expected discoveries. 

Mallen-Ornelas et al. (2001) expected to find 1 
transiting (V)HJ. Although this is very low number 
statistics, this value is in good agreement with our result 
for Pi,is2 — 0.1. The absence of any detection is also 
compatible with our estimates. We point out that good 
weather conditions could have led to the discovery of a 
few {V)HJ, showing that searches using 4m-class tele- 
scopes and a wide-field camera are promising strategies. 
EXPLORE-I lasted only 11 nights. A simple scaling of 
the results for a 3-months survey shows that a harvest of 
several tens of {V)HJ is not unrealistic. 



Weather conditions 


VHJ 


HJ 


2 0.5 


Real conditions (fc = 2) 


0.4 (0.5) 


1.2 (1.5) 




Good weather (k > 2) 


3.4 (4.2) 


2.8 (3.4) 





Good weather (A; > 3) 


2.2 (2.7) 


0.7 (0.9) 
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Fig. 5. Distribution of the expected number of {V)HJ for 
EXPLORE-I (with good weather conditions) as a function 
of the number of transits observed {upper left), of the dis- 
tance {upper right), the magnitude {bottom left) and the 
spectral subtype {bottom right) of the host star. The ir- 
regular shape of the histogram as a function of distance is 
an artificial effect caused by the rounding off of the stellar 
magnitudes to the nearest integer or half-integer. 

3.3. OGLE-III 

The Optical Gravitational Lensing Experiment entered its 
third phase, OGLE-III, in June 2001. It took place at the 
Las Campanas Observatory, Chile, using the 1.3m Warsaw 
telescope and the 8k MOSAIC camera, with a total field 
of view of 0.34 deg^. All observations were made through 
the 7 filter. Four surveys have been carried out so far : 

- OGLE-III-1 (June 12 to July 28, 2001). More than 
800 images of three fields in the direction of the galac- 
tic bulge were collected within 32 nights. The exposure 
time was 120 seconds, and each field was observed ev- 
ery 12 minutes. 

- OGLE-III-2 (February 17 to May 22, 2002). More 
than 1100 images of three fields located in the Carina 
region of the galactic disk were collected in 76 nights. 
The exposure time was 180 seconds, and the temporal 
resolution was about 15 minutes. 

- OGLE-III-3 (February 12 to March 26, 2003). The 
photometric data were collected during 39 nights span- 
ning the 43 days of the survey. Three fields of the galac- 
tic disk were observed with a time resolution of about 
15 minutes. The exposure time was 180 seconds. 

- OGLE-III-4. Starting on March 25, 2003, this survey 
collected its main photometric material until middle 
May 2003, but observations were also collected until 
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Fig. 6. Distribution of the {V)HJ discoveries for OGLE- 
III-2 as a function of the number of observed transits {up- 
per left), distance [upper right), magnitude {bottom left) 
and of the host star spectral subtype (bottom right). 



July 25, 2003 with sparse temporal sampling. Three 
fields of the galactic disk were observed with the same 
exposure time and resolution as the two previous sur- 
veys. 

Using the observation windows and instrumental pa- 
rameters of these 4 surveys (A. Udalski, private commu- 
nication), we have computed the expected harvests for 
the two detection criteria (see Table OJ. OGLE-III-1 and 
OGLE-III-2 have yield 5 genuine extrasolar planets so far. 
OGLE-III-3 and OGLE-III-4 have discovered 40 transit- 
ing companions, but follow-up spectroscopy has not given 
the final harvest yet. We note that, although predicting 
absolute planet counts is not our main goal, our OGLE- 
III expected harvests are in very good agreement with the 
actual ones. Our work hypotheses thus allow to compute 
realistic harvests in the case of {V)HJ. 

Fig.Elshows the distribution of the {V)HJ discoveries 
predicted by our simulation for OGLE-III-2 as a function 
of the characteristics of the host star and of the number of 
observed transits. Most detections occur for stars located 
between 1500 and 2500 pc from the Sun, well matching 
the range of distances for the real planets OGLE-TR- 
111 , OGLE-TR-113 and OGLE-TR-132, i.e., about 1500 
pc ( [Schneider, Extrasolar Planet GataIog| ). We also note 
that our predicted ratio of VHJ to HJ is in good agree- 
ment with the observed one (see Table Future planet 
harvests with OGLE-III-3 and OGLE-III-4 shall allow 
to better constraint the absolute counts of VHJ around 
LMSS. 



Table 3. Results predicted by our simulations for OGLE- 
III surveys, compared to the actual harvests. The values 
are given for f3 = 9 and the values in parenthesis cor- 
respond to /3 — 7. No actual harvest is available yetfor 
OGLE-III-3 and -4 yet. 



Survey 


Planet 


Predicted 


Actual harvest 


OGLE-III-1 


VHJ 


1.6 (2.2) 


1 




HJ 


0.7 (0.9) 


1 


OGLE-III-2 


VHJ 


2.7 (3.6) 


2 




HJ 


1.5 (2.0) 


1 


OGLE-III-3 


VHJ 


2.1 (2.8) 






HJ 


0.8 (1.1) 




OGLE-III-4 


VHJ 


2.6 (3.6) 






HJ 


1.8 (2.5) 





3.4. Analysis of fictitious surveys 

3.4.1. Setting the parameters of the surveys 

All the surveys we consider in the following use red filters, 

1. e., redder than the y-band. This choice is motivated by 
the fact that limb darkening and atmospheric and galac- 
tic absorptions are minimized in the red. It also maxi- 
mizes the number of available LMSS, hence increasing 
the probability of detecting extrasolar planets. Extending 
the above argument from the visible to the near-infrared, 
we have included surveys using the J, H and K filters, 
on VISTA-IR. However, we have not included the longer 
wavelength filters such as L or M, because the quantum 
efficiency of the detectors is much lower with present-day 
instrumentation. 

The choice of the exposure time is a critical issue for 
surveys carried out in the visible. Long exposure times 
lead to high SNR, but also to a far too large number of 
saturated stars in galactic fields. Increasing the exposure 
times also increases crowding, and consequently decreases 
the SNR for any contaminated star. As mentioned in Sec. 

2, effects of crowding are not taken into account here, and 
will be the subject of a forthcoming paper. If the goal 
of the survey is to discover planets as small as possible, 
a compromise has to be found between the duration of 
the exposure time plus readout time of the CCD, the de- 
sired SNR in each individual exposure and the number 
of sampling points during an eclipse. We have estimated 
the exposure time leading to the largest number of unsatu- 
rated LMSS with sufficient SNR and chosen to test three 
different exposure times close to this "optimal" value for 
each optical instrument. The actual readout time of each 
instrument is also taken into account in our calculations. 

The situation is simpler for near-IR instruments, since 
the readout time is negligible. The harvest is thus largely 
insensitive to the adopted exposure time and we have 
adopted the exposure time leading to the largest number 
of unsaturated LMSS. 

Finally, we have tested surveys of four different dura- 
tions: 30, 60, 120, and 180 nights. 
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Fig. 7. Distribution of predicted telluric planet discoveries 
(/? — 7) for VISTA-IR and the 4-years J-band survey, as 
a function of the number of observed transits [upper left), 
distance [upper right), magnitude [bottom left) and of the 
host star spectral subtype [bottom right). 

We have checked that there exist fields in the galactic 
disk that can be observed continuously during six months, 
either from Paranal or Mauna Kea. However, their visibil- 
ity is reduced to about 4 hours during the first and the 
last month of a six months period. We have used the sea- 
sonal observability of typical fields of the galactic disk to 
build our window function and we have also investigated 
the interest of extended surveys, carried out over several 
years (between 1 and 4 years) , similarly to space missions 
such as COROT and KEPLER, but with non-continuous 
visibility. 

The predicted harvests for the different surveys are 
summarized in Table 21 The columns in all Tables are 
labeled as exposed in Section 2. 

3.4.2. Results 

In the following, we comment on the results presented in 
Table H 

— Optical filters: optical surveys are about twice as 
efficient in the /-band as in the i?-band. This is mainly 
due to the increased brightness of most LMSS in the 
/ filter, which more than compensates for the lower 
throughput of the instruments in the /-band. For this 
reason, we consider only /-band searches for surveys 
longer than I month. 

— Near-IR filters: infrared surveys are more efficient 
in the J- and //-bands than in the /sT-band. In the 



Fig. 8. Same as Fig. but for VISTA- Vis and 4 years of 
observation through the / filter. 

infrared, the gain in the number of LMSS observed 
in K does not compensate for the lower sensitivity of 
the detectors at this wavelength. The J filter appears 
to be the best choice. 

— Exposure times: in the optical, the longer exposure 
times do not always yield to larger harvests. Although 
this increases the SNR across the duration of the tran- 
sit, it also leads to more saturated stars, so that the 
improvement in the harvest is negligible. 

— (V)HJ vs. IZ planets: even for the shortest survey 
duration (30 days) and for the least favorable detection 
criterion [(3 = 9), the harvests in [V)HJ are always 
large, with tens of discoveries. Increasing the duration 
of the survey does not increase very much the number 
of VHJ detected, but increases drastically the number 
of giant planets found in the HZ and in the IZ , the 
largest gain being for surveys longer than one year. 
This is simply due to the better matching of the time 
base-line with the longer revolution periods of planets 
in the HZ and the IZ. 

— Field vs. depth in the optical: among the exist- 
ing optical instruments, the harvests obtained with 
SUBARU are smaller than those obtained at the 
CFHT, despite the larger telescope size: the improved 
depth of SUBARU does not compensate for the smaller 
field of view. 

— Telluric planets: Under the assumptions used here, 
our fictitious ground-based surveys seem to have a 
low potential for discovering TIZ . VISTA- Vis has the 
highest potential for this purpose, followed by VISTA- 
IR. Habitable planets seem out of the reach of ground- 
based surveys. The telluric planets which could be dis- 
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Table 4. Results obtained for fictitious surveys. ET = exposure + readout time. The values are for /3 = 9. The values 
given between parenthesis correspond to /3 = 7. 
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Fig. 9. Distribution of predicted planets (all types in- 
cluded) discoveries using /3 = 7 for the COROT mission, 
as a function of the number of observed transits (upper 
left), distance [upper right), magnitude [bottom left) and 
of the host spectral subtype [bottom right). 



covered by VISTA-Vis or VISTA-IR would probably 
orbit around M-dwarfs relatively close to the Sun (see 
Figs. Hand El). 



3.5. Expected harvest of the COROT mission 

COROT (convection, Rotation and planetary Transit) 
will be the first satellite launched with the aim of detecting 
exoplanets using the transit method l |Rouan et al., 2000| ) . 
Another goal of this mission is to carry out several projects 
in asteroseismology ( |Baglin et al., 1998| ). It uses a 27 cm 
telescope in combination with two 2048x2048 CCD cam- 
eras, and a final field of view of 3.5 deg^. During the 2.5 
years of the mission, 5 fields in the galactic plane will be 
observed continuously, each one during 150 days. The ex- 
oplanet detection capability of this mission has already 
been analysed ( |Borde et al., 2003| ) but, since our goal is a 
comparison between surveys, we analyse the planet har- 
vest expected from COROT using the same assumptions 
as for the other surveys considered. 

COROT will observe fields at low galactic latitudes, 
i.e., in directions with a high density of LMSS. Indeed, 
each field will contain up to 12,000 dwarf stars with visual 
magnitudes between y = 11 and V = 16.5. The exposure 
time will be 31.7 seconds, without any filter. 

We use the expected number of photoelectrons per 
exposure for every magnitude bin between V = 11 and 
= 16.5 as well as the different noise contributions given 




Fig. 10. Magnitude distribution of the V < 16.5 LMSS 
in a typical COROT field. 

Table 5. Results obtained for COROT, integrated over 
the duration of the whole mission. 



/3 


VHJ 


HJ 


GIZ 


TIZ 


9 


6.46 


13.94 


1.78 


0.01 


7 


7.59 


17.24 


2.27 
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Table 6. Results obtained for the KEPLER mission. 
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on the COROT web site'^ and in Borde (2003) to compute 
the theoretical SNR per exposure. Using in addition the 
duration of one field observation (150 days continuously), 
and multiplying the obtained harvest by 5 to take into 
account the 5 fields that will be monitored, we compute 
the harvests presented in Table|5| The distributions of the 
total number of planets detected for /3 = 7 as a function 
of the characteristics of the host stars and of the number 
of transits observed are shown in Fig. [HI Fig.llOlshows the 
distribution of the LM SS as afunction of their magnitude 
in the F-band, for a typical field observed by COROT. 

Borde et al. (2003) estimated that COROT wiU be 
able to detect numerous giant planets, a prediction that 
is confirmed by our results. It is nonetheless immediately 
apparent that the expected number of giant planets dis- 
coveries (~ 25) after 2.5 years is well within the reach of 
30 days ground-based surveys. Only SUBARU would need 
a survey of several months to be able to compete, due to 
its smaller field of view. 



http: / / www.astrsp-mrs.fr /projets /exoplan/corot / 
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Fig. 11. Distribution of predicted giant planet discoveries 
using /3 = 7 for KEPLER, as a function of the number 
of observed transits {upper left), distance {upper right), 
magnitude {bottom left) and host star spectral subtype 
{bottom right) . 



Fig. 13. Same as Fig.HIlbut for THZ (/3 = 7). 

ble, even though we only consider telluric planets similar 
to the ones in our Solar System. 
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Fig. 12. Same as Fig. E] but for TIZ (/? 7). 

Borde et al. (2003) argue that COROT wiU be able 
to detect large terrestrial planets {R > i?0, with i?^ — 
earth radius) in close orbits, provided that this kind of 
object is frequent. Under our assumptions, however, the 
probability that COROT finds a telluric planet is negligi- 



3.6. Expected harvest of the KEPLER mission 

The KEPLER mission is dedicated to the detection 
of earth-like planets orbiting sun-like stars in the HZ 
dKoch et al., 1 998). This satellite is scheduled for launch 
in 2007, and will observe continuously one single 105 deg^ 
field during 4 years. It will monitor ~ 100, 000 main se- 
quence A-K stars with a mean magnitude of y = 14 using 
a 0.95 m telescope and 3 seconds exposures. Using the 
technical data available on the KEPLER web site^, we 
compute the SNR for an individual exposure time, and 
use it to estimate the planet harvest of KEPLER under the 
same assumptions as in the previous sections. KEPLER 
will observe a field with a mean galactic latitude of 13.3°, 
leading to a maximal distance Dl ~ 1300 pc according to 
Eq. n The results are presented in Table IHl Figures E] 
[T7\ and ^1 show the distribution of expected planet dis- 
coveries {p = 7) as a function of the number of observed 
transits and of the characteristics of their host star. 

KEPLER should be able to detect many telluric plan- 
ets, even in the HZ. In contrast, the expected harvest in 
giant planets is smaller than the one expected for 30-days 
surveys carried out with VISTA- Vis. This is due to the 
larger field of view of KEPLER and to its shorter expo- 
sure times. This strategy is certainly the best one for hunt- 
ing telluric planets, but does not allow to observe enough 
LMSS to detect a large number of giant planets in close 
orbits. 

littp://www.kepler. arc. nasa.gov 
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Total = 35298 
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V 

Fig. 14. Magnitude distribution of the V < U LMMS in 
the KEPLER field of view. 

According to the KEPLER web site, the expected re- 
sults of the mission are the following: 

— assuming that most telluric planets have R ~ 1-R©, 50 
discoveries should be made. This number increases to 
185 if most telluric planets have R ~ 1.3i?© and to 640 
if most of them have R ~ 2.2Rq). 

— 135 inner-orbit giant planets and 30 outer-orbits giant 
planets should be discovered by KEPLER, plus an- 
other 870 (y)HJ discovered through the modulation 
of their reflected light. 

These expectations are based on the following assump- 
tions : 

— The number of main-sequence stars monitored is ^ 
100, 000, after the exclusion of the most active dwarf 
stars. 

— The variability of these ~ 100, 000 stars on the time- 
scale of a transit is solar or close to solar. 

— All stars host in average two earth-size planets in the 
region between 0.5 and 1.5 AU. 

— Every star host one single giant planet in a jovian-like 
orbit. 

— On average, 1% of the target stars have a giant planet 
with an orbital period smaller than 1 week, 1% with 
periods between 1 week and 1 month and 1% with a 
period between 1 month and 1 year. 

Detections around binaries are taken into account in 
these expectations. The detection criterion used by the 
KEPLER team is A^„i„ = 3, = 8. 

Our estimate of the number of LM SS in the KEPLER 
field leads to a value of ~ 35,000. Fig. 1141 shows the dis- 
tribution of these LMSS as a function of their apparent 
magnitude. The discrepancy between our number of target 
stars and the one from the KEPLER team may have two 



origins: (1) we do not take into account the spectral types 

0. B and A, (2) we assume a constant density of LMSS 
in a given volume, and nothing outside, (3) the actual 
extinction coefficient Ay of the KEPLER field should be 
somewhat lower than 0.7 mag/kpc. 

For the giant planets, scaling our result (the sum of 
the first 4 columns in Table (HJ ~ 33 planets) to the num- 
ber of targets stars expected by KEPLER and taking into 
account the possible detection around binaries would lead 
to a value in good agreement with the one of the KEPLER 
team. 

Our assumptions about telluric planets are quite dif- 
ferent from the ones used by the KEPLER team, but our 
results 50 — 80) are in good agreement, provided that 
the values presented on the KEPLER web site include all 
telluric planets, and not only habitable ones. 

KEPLER is a key project in the search for life outside 
our Solar System. We can remark from Fig. El that hab- 
itable planets would be detected mainly around G and K 
dwarfs located about 200 pc from the Sun, thus around 
nearby solar-type stars. The predicted number of THZ 
(2 — 5) is small. It could become much higher if our radii 
distribution proved to be too pessimistic. 

4. Discussion - Conclusions 

The main purpose of our simulations was to weight the ad- 
vantages and disadvantages of various ground- and space- 
based searches for exoplanets using the transit method. 
Tables 2 to 6 summarize our results and provide the ex- 
pected harvest in exoplanets for a broad variety of tele- 
scope/instrument combinations. Our main conclusions are 
the following: 

1. As far as telluric planets in the Habitable Zone are 
concerned, space-based surveys are the only viable op- 
tion. Such searches remain extremely difficult. They 
not only require a space instrument, but also a very 
wide field of view. From space, only the KEPLER 
mission should be able to find telluric planets in the 
Habitable Zone. 

2. Telluric planets in the Intermediate Zone are much eas- 
ier to discover. KEPLER could detect more than 40 
of them during its four years of observations. Only a 
few (1-2) TIZ might be discovered from ground-based 
surveys of the same duration, using VISTA- Vis. 

3. Ground-based searches are better than space searches 
at finding giant planets. While KEPLER is about 
as efficient as CFHT at finding giant planets in the 
Habitable Zone (with 1 expected discovery vs 2), a 
CFHT search is 4 times better than KEPLER at find- 
ing the same planets in the Intermediate Zone, and 
5 times better for {V)HJ. This is due to the much 
deeper exposures. 

Ground-based and space-based transit searches are 
complementary. Because they go deeper, ground-based 
searches easily find large planets with a short period, 
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such as (V)HJ. Space searches remain mandatory for 
telluric planets. COROT will be a pionner mission in 
this field, but we shall have to wait for KEPLER to 
obtain a significant harvest of such objects, provided 
that they are not significantly less common than expected. 

The above results give orders of magnitude estimates 
for the expected harvests and allow to emphasize the 
relative merits and drawbacks of the different searches. 
However, a word of caution should be given, to avoid over- 
interpretation of the results. One need to be aware that: 

— Space missions like COROT or KEPLER will defocus 
the images, increasing drastically the size of the Point 
Spread Function (PSF). While this will minimize sat- 
uration of bright objects and increase the SNR per 
image, it might result in severe image blending. The 
corresponding loss of efficiency in transit detection will 
depend on the method used to post-process the data, 
as do the effects of blends. This is why we have deliber- 
ately chosen not to take PSF convolution into consider- 
ation and to leave it for a future work. All our estimates 
are therefore upper limits on the expected harvests. 

— The weather simulations used are very simple, and 
could be somewhat optimistic. Real weather conditions 
could lead to lower harvests for ground — based sur- 
veys. 

— The photometric techniques required for the analysis of 
ground-based surveys are very efficient in the optical 
(e.g., PSF fitting, image subtraction, image deconvolu- 
tion) , but may be less efficient for near-IR data, where 
the sky subtraction is more critical. We therefore ex- 
pect our near-IR harvest estimates to be slightly more 
optimistic than the optical estimates. 

— In the near-IR, second-order extinction effects can 
have a large impact on the photometric accuracy 
( [Bailer- Jones fc Lamm, 2003| ). Time dependent atmo- 
spheric extinction depends on the spectral energy dis- 
tribution of the target. This may imply the need to use 
reference objects of the same spectral type as the tar- 
get stars and complicate further the analysis of near- 
IR data. 
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